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I calculate the self-force acting on a particle with electric charge q moving on a generic geodesic 
around a Schwarzschild black hole. Using methods similar to those developed for the scalar field 
case discussed in fl|], 1 investigate the relative sizes of the conservative (half-advanced plus half- 
retarded) and dissipative (half-advanced minus half-retarded) pieces of the self-force. I also display 
the regularization parameters used in the mode-sum regularization scheme. 
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I. INTRODUCTION 

This is the second paper of a series of papers study- 
ing the self-force on a point particle in generic geodesic 
orbit around a Schwarzschild black hole. I extend the 
previous calculation of the scalar self-force [S] to elec- 
tromagnetism, studying in particular the effects of the 
conservative part of the self-force. 

A test particle in orbit around a black hole will follow 
a geodesic. Going beyond the test mass limit, this is no 
longer true and the particle's path will deviate from a 
geodesic of the background spacetime. As seen from the 
background spacetime, the particle is said to experiences 
a self-force due to its interaction with its own field. In or- 
der to accurately model the motion of the body, including 
its inspiral toward the black hole, I seek to evaluate the 
self-force and calculate its effect on the motion. Several 
methods to achieve this have been proposed in the liter- 
ature I elect to use the mode-sum regularization 
scheme introduced by Barack and Ori [3|], which been 
proven to be highly accurate. 

In this paper, rather than dealing with the gravita- 
tional problem, I focus on the technically simpler prob- 
lem of a point particle endowed with an electric charge 
q orbiting a Schwarzschild black hole of mass M . In this 
context I use a numerical simulation to check the ana- 
lytically calculated regularization parameters used in the 
mode-sum regularization scheme, which I calculate in a 
manner analogous to 0. This calculation also makes it 
possible to investigate the behaviour of the conservative 
(half-advanced plus half-retarded) part of the self-force in 
the strong- field limit, extending previous work by Pound 
and Poisson Q. Different from the scalar case calcu- 
lation, where the conservative self-force is suppressed, 
the conservative electromagnetic self-force appears at the 
same post Newtonian order as the gravitational conser- 
vative self-force. Agreement, even if only qualitative, be- 
tween the results for the electromagnetic problem, where 
our physical intuition allows us to understand the mech- 
anisms at work, and those for a point mass recently ex- 
plored by 0, Q can thus help provide a clearer under- 
standing of the mechanisms at work in the gravitational 



case as well. 

Throughout the paper I use geometrized units in which 
G = c = 1 and the sign conventions of ^] . 



A. The problem 

Since my approach is essentially identical to that de- 
scribed in and [l| (paper I and paper II from now on) , 
I will only briefly introduce the required notation. 

The first order self-force is calculated on a geodesic 
of Schwarzschild spacetime, whose metric is written in 
Schwarzschild coordinates as 



ds^ = -/di^ + /"Mr^ + r'^dVL, 



(1.1) 



where / = (l-^), dVl = [dO^ + sit? eAcjP) is the 
metric on a two-sphere, and t, r, 9 and are the 
usual Schwarzschild coordinates. I numerically solve the 
Maxwell equations 



g^^^^Fapix) = 47rj„(a;), 

ja{x) Ua{T)S4{x,z{T))dT, 
J J 



(1.2) 
(1.3) 

(1.4) 



where Vq is the covariant derivative compatible with the 
metric (7a/3, Fap is the Faraday field tensor sourced by 
a charge q which moves along a world line 7 : r i— > 
z(t) parametrized by proper time r. The current den- 
sity jaix) appearing on the right-hand side is writ- 
ten in terms of a scalarized four-dimensional Dirac S- 
function ^4(x, x') = S{xo — x'q)S{xi — x'i)S{x2 — x'2)S{x3 — 
Xg)/-^/— det(gQ,^). After having obtained the Faraday 
tensor I regularize it using the mode-sum regularization 
scheme introduced by Barack and Ori [Sj 
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MM 



(1.5) 
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where indices in parenthesis (p) signify components with 
respect to an orthonormal tetrad e" ^) and the coefficients 

^^'^ ^"^^ independent of 



£; they are hsted in Appendix [Bj Finally I compute the 
regularized self-force 



(1.6) 



from the regularized Faraday tensor and the four velocity 
of the particle. 



B. Organization of this paper 

In Sec. [TT] I introduce the ideas behind the discretiza- 
tion scheme used in the numerical simulation. Sec. IIIII 
describes the choices I make in order to handle the prob- 
lems of specifying initial data and proper boundary con- 
ditions. In Sec. I VIII I describe the tests I performed in 
order to validate my implementation of the numerical 
method. Sec. I Villi contains sample results for a small 
number of representative simulations. Finally in Sec. IIXI 
I calculate the conservative self-force for the same set 
of simulations. The appendices contain technical details 
and an alternative calculation using the vector potential 
instead of the Faraday tensor. 



II. NUMERICAL METHOD 



In this section I describe the algorithm used to 
integrate the Maxwell equations numerically. I use 
the second-order algorithm introduced by Lousto and 
Price ^0|] suitably extended to handle a coupled system 
of equations. 



A. Wave equations for the Faraday tensor 

I introduce a vector potential Aa in terms of which the 
Faraday tensor is given by 



F 



a/3 



A 



f},a 



(2.1) 



where a comma denotes an ordinary derivative. I use 
vector spherical harmonics = OaY^"^ and — 
ca^ObY^"^ , where cab is the Levi-Civita tensor associ- 
ated with the metric ^Iab on the two-sphere (eg^ = sin0), 
as introduced in [ll|, [l3| . I decompose the vector poten- 



tial and the current density into 

v4<,(^,r,^^,0) = A^(^,r)r,™(0,</.), 

AA{t,r,e,c^) = vtUt.r)Z'r{e.cj)) 
jAit,r,0,<p)^jr^'^{t,r)Z'/^{9,<l,) 



(2.2a) 

for a = t,r, 

(2.2b) 

(2.2c) 



for A = 0, 



(2. 2d) 



where a summation over £ and m is implied. Substitut- 
ing these into Eq. (|1.2I) I arrive at two sets of coupled 
equations for the even (A^™, v^rn) and odd (vim) modes 

J o 9 ^ ^ 



Qj,2 



dtdr 



dr 



dt 



^a^^^ (2.3a) 
at 



dt'^ dtdr 



dr 



(2.3b) 



,dA: 
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Qj.2 

2M 



dr 



dt 
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even 
m ' 



2M dv. 



em 



+ 



dt^ dr^ 
+ 1).. 



J.2 Qj, 



47rj, 



odd 

e-m ' 



(2.3c) 



(2.3d) 



where 



■em 



-^^^"(^,^0)5(^^0), (2.4a) 



fr = ^Y^"'i^,^omr-r,), 



■even 

Jem 



■odd 

Jem 



imqfj -^^^ 



(2.4b) 



Y'"^{^,^o)Sir-ro), (2.4c) 
.lLL_d,Y'"\^,^o)S{r~ro). (2.4d) 



£{£+l)Er^ 



In the equation above an overbar denotes complex con- 
jugation, an overdot denotes differentiation with respect 
to r, £^ = —Ut is the particle's conserved energy per unit 
mass, J — its conserved angular momentum per unit 
mass, and u" — is its four velocity. Quantities bear- 
ing a subscript "0" are evaluated at the particle's posi- 
tion; they are functions of r that are obtained by solving 
the geodesic equation 



u'^Vpu"' = 



(2.5) 



in the background spacetime. Without loss of generality, 
I have confined the motion of the particle to the equato- 
rial plane = f . 
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The three even mode equations Eq. (I2.3a|) - Eq. (|2.3cp 

are not yet amenable to a numerical treatment, as they 
are highly coupled. In order to obtain a more convenient 
set of equation I define the auxiliary fields 







\ dr 


dt 




dr J 




dt ' 



(2.6) 
(2.7) 
(2.8) 

which, up to scaling factors, are just the even multipole 
moments of the tr, r<j) and tcj) components of the Faraday 
tensor 

FtA = ^(-C'" + if XT). (2.10) 

(2.11) 

e.m 

= -^^(^ +!){)£„ sin(0)r^". (2.12) 
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I note that the three fields i/)^™, x^™ and are not 
independent of each other, in fact knowledge of -0^™ is 
sufficient to reconstruct x^™ and Eq. (|2.3ap can be 
rearranged to yield 



f df 



im 



An 



£{£+!) dr £{£+!) 
and similarly from Eq. ()2.3b|) 



Jt 



Em 



(2.13) 



1 



+ dt + 



■em 

Jr ' 



(2.14) 



showing that knowledge of ■0^™ is sufficient to reconstruct 
the even multipole components of the Faraday tensor. In 
this work however I choose to solve for and di- 
rectly, rather than to numerically differentiate -0^™ to ob- 
tain them. The gain in speed from reducing the number 
of equations does not seem to offset the additional time 
required to calculate tp^^ accurately enough to obtain 
good approximations for its derivatives at the location of 
the particle. In this approach Eqs. (|2.13p and (|2.14p are 
treated as constraints that the dynamical variables have 
to satisfy. 

Dropping the superscripts i, m for notational conve- 
nience and following fl^ I form linear combinations of 
derivatives of Eqs.(I23a| - ([^^c)) . I use [dr{r^ (EjEJ) - 



(9t(r2 (HJil))] for and find 



^0 - 47r/ 



dr 



dt 



(2.15a) 
(2.15b) 



where V = l{t + 1)'-^^ and r* = r + 2M\n{^ - 1) 
is the Regge- Wheeler tortoise coordinate. Similarly I use 
[/ (l2:3bl) - dr{f (I23cl) )] for x and [ (l23al) - dt^M] for 
t I find 



Qj,*2 Q-f.2 



-Vx = Sx, 

ffr" 



dr 



d^£ d^£ 



dt 



(2.15c) 
(2.15d) 

(2.15e) 
(2.15f) 



where = 2{r -zM){r -2M) ^ ^j^jjg g^j^ partially coupled 
Eqs. (j2.15bl) - (j2.15f|) are much easier to deal with than 
the original set Eqs. (I2.3ap - (|2.3cp . The coupling is in 
the form of a staggering, which allows me to first solve 
for -0 and use this result in the calculation of ^. On the 
other hand, the source terms appearing on the right-hand 
side contain derivatives of Dirac's (5-function resulting in 
fields that are discontinuous at the location of the par- 
ticle. Lousto's scheme is designed to cope with precisely 
this situation. 

I derive explicit expressions for the source terms Sa on 
the right hand sides 



47rq /.. imrQj\- tt Motrin 
-^fo [ro -3— ) y£„(-,(po), (2.16b) 



Sa = Ga{t)foSir - ro) + Fa{t)f5'{r - r„), (2.16a) 
G^it) 

F0(i) =47rg/o ' ''° 
Fxit) = 



.F^ 
47rqfo 



2 ^]Ytm{-,y^o) 



TT 



Aug Jim n 
"i?^(f+l)rg^o^''"^2'^°^ 



(2.16c) 
(2.16d) 
(2.16e) 



G^{t) = -Anq 



Jim 



'/2M 


2fo\ 




ro J 



E^i{e + l)r2 

2 f'/o^^m(-;r,<Po), 



Fdt) 



AnqJimfo 2t 

F^e{e + i)r2 ^° "'"^2^ 



(2.16f) 



My functions Ga and Fa correspond to G/fo and F/f 
in [l3|, respectively, they are independent of r (but do 
contain terms in ra(t)). I prefer this form of the source 
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terms over the form given in jl^ since it simplifies the 
integral over the source term Eq. (3.6) of (loj 



dAS 



G(ro{t),t) 
l-2M/roit) 
F{r,t) 



dr\l- 2M/7 



r=ro(i) 



dt 



^^[l-2Af/ro(<2)]2^ ■ ^ ' 

Since 0^°"^*° = /oG„(<) and 6^°""*° = the 
first term in square brackets inside the integral simpli- 
fies, while the second term vanishes completely. Fa only 
appears in the boundary terms. 



B. Constraint equations 

The full set of Maxwell equations consists of the in- 
homogeneous equations Eq. (jl.2p as well as the homo- 
geneous constraints Eq. (|1.3p which have to be satisfied 
by a solution to Eq. (|1.2p . In the usual approach intro- 
ducing a vector potential implies that the constraints 
are identically satisfied since they reduce to the Bianchi 
identities for the second derivatives of Aa ■ When solving 
for the components of the Faraday tensor directly there 
is no a priory guarantee that a solution to Eq. (|2.15b|) - 
(|2.15f|l . and (|2T3d|) satisfies Eq. It turns out, how- 

ever, that a decomposition into spherical harmonics is 
sufficient to show that all but one of the constraints are 
identically satisfied. The one that is not identically true 
is the trip (or trO) equation, which in terms of -0, X and 
^ reads 



,2 f +4.' 



(2.18) 



If the fields satisfy the sourced Maxwell equations 
Eqs. ((23a|) . (|2.3bp . then Eq. (|2T8| is just the evolution 
equation for -0. Thus Eq. (|2.18p is valid whenever il) sat- 
isfies the consistency relations Eq. (|2.13p and (|2.14l) . 

Analytically then, the situation is clear. Given a set of 
compatible initial conditions for i/), x and ^ which initially 
satisfy the constraint equations, a solution to the system 
of Eq. (|2.15bp - ((2T5^ . ((23d)) satisfies the full set of 
Maxwell equations at all later times, too. 

Numerically I monitor but do not enforce Eq. (|2.13p 
and (|2.14l) . I generally find that violations of the con- 
straints are at least three orders of magnitude smaller 
than the field quantities themselves. Figures [1] and [5] 
compare x obtained from its evolution equations to that 
obtained from Eq. (|2.14p . 
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FIG. 1: Violations of the constraint = 
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A I e(e+i) at 

in the vacuum region away from the location of particle 



= 
I 



plot the X ™d logio l^xl ^ obtained on a spatial slice at 
time t = 600 M. For this slightly eccentric orbit (p = 7.0, 
e = 0.3) using a stepsize h = 1/512A/ the errors in the 1 = 2, 
m = 2 mode are at least three orders of magnitude smaller 
than the field values. The exponentially growing signal be- 
tween 300A/ < r* < 500 is a remnant of the initial data pulse 
travelling outward. 
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FIG. 2: Violations of the constraint 







— X + 1(1+1) at 

at the location of the particle as a function of time. I display 
X and logjQ \Z-x\ for the ^ = 5, m = 3 mode of a particle 
on an eccentric orbit with p = 7.8001, e = 0.9 with stepsize 
h = 1/256 M. During the time 400 M < t < 800 M the 
particle is in the whirl phase. The exponentially decaying 
signal before t ~ 250M is the initial data pulse. 



C. Monopole mode 



For the electromagnetic field, the monopole mode t — 
is non-radiative. The vector harmonics Z^J^ and X^'" 
cannot be defined in this case and the only surviving 
multipole mode is ■0. For the monopole case Eq. (l2T5b)) 
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reduces to a wave equation in flat space 



dr* 



dr 



dt 



(2.19) 



which is simple enough so that I can solve it analytically. 
A straightforward calculation shows that 



(2.20) 



satisfies Eq. (|2.19p and corresponds to no outgoing radi- 
ation {dt — dr*)ip = at the event horizon and no ingoing 
radiation (dt + dr*)^ = at spatial infinity. 



D. Discretization — even sector 
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Lousto's method is directly applicable to terms of the 
form — + ^-^i V{r)^l: (ie. the wave operator and po- 
tential terms) on the left-hand side of the equation and 
the source terms Sa{t) on the right hand side. Here i}} is 
used as a placeholder for any one of '0, x or ^; V{r) is an 
expression depending only on r. I discretize these as 



^"^^ 1 - + ) = + ^2 - 01 - 04] , 

(2.21) 

0{h'^) vacuum cells 
0{h^) sourced cells, 
(2.22) 



du dv V{r)ip 



dr* 



and 



du du Sa{t) 
cell 



= 2 



Gc,{t,ro{t))dt 



± 



l-2M/r(ti, 
2F^{t2,ro{t2)) 
1 - 2M/r{t2) 



[l±fo{t2)/EY 



(2.23) 



where u = t—r*, v = t+r* are null coordinates, ipii- ■ ■ :0'4 
refer to values of the field at the points labelled 1,. . . ,4 in 
Fig. ^ h = At = Ar* /2 is the step size, Vq is the value of 
the potential at the centre of the cell, Ai,. . . ,^4 are the 
areas indicated in Fig. |3] and ti and ^2 are the times at 
which the particle enters and leaves the cell, respectively. 

Spelled out explicitly the evolution equations for vac- 
uum cells are 
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+ (1- Y^o)(V'iH 




(2.24a) 


X3 = 


-X2 


+ (1 - Y^o)(xi ^ 
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(2.24b) 
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(2.24c) 



FIG. 3: Points used to calculate the integral over the potential 
terms. Grid points are indicated by blue circles. 



and for sourced cells 

"03 



[1 + -j(^2 - + [1 - -^(^4 + A3)]V4 



[1- -j(^l+A3)]0'l 

i(l-^A3) [(dudvS^it), 



(2.25a) 



X3 = -[1 + ^(^2 - A3)]X2 + [1 - -j(A4 + ^3)]X4 



[1- Y(^1+^3)]X1 
1(1-^^3) ffdudvS^it), 



(2.25b) 



6 = -[1 + ^{A2- A3)]6 + [1 - -^(^4 + ^3)1^4 

+ [i--j(^i + ^3)]a 

- iv5,o(Al0'l + ^2V'2 + ^3V'3 + ^40^4) 
-1(1-^^3) ffdudvS^it). 



E. Discretization— odd sector 



(2.25c) 



When written in terms of r*, Eq. (j2.3dp . which governs 
the odd modes ■0^™, is 



dHim _ d^Vim _ £{e+l){r- 2M) ^ 
dr*^ dt'^ r^ 



31 



dd 



(2.26) 

^^-^deY'^C-^.^^)8r -rl). (2.27) 



l(i + \)Erl 



Eq. (|2.26p is of the form of the scalar wave equation 
discussed in paper II. I re-use the fourth order numer- 
ical code described there with V = ^(^+i)(^~2^^) ^ g — 
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An 



qfJ 



:^doY^^{^,ipo). This yields accurate results 



for V and its derivatives. 



values for ip and dr*ip. Specifically I fit a fourth order 
polynomial 



III. INITIAL VALUES AND BOUNDARY 
CONDITIONS 

I follow the approach detailed in paper II for the scalar 
self-force and do not specify physical initial data or an 
outer boundary condition. I arbitrarily choose the fields 
to vanish on the characteristic slices u = uq = to — Tq 

and V — Vq = Iq + rX 



ip{u = uq) = ^{v = Vq) = 0, 



(3.1) 



thereby adding a certain amount of spurious waves to the 
solution which show up as an initial burst. 

I implement ingoing wave boundary conditions near 
the event horizon, sufficiently close that numerically 
r 2AI, so that the potential terms in Eqs. (j2.15bp - 
(|2.15fl) vanish. This happens at r* w —73 A/ and I im- 
plement the ingoing waves condition dutp — there. Near 
the outer boundary this is not possible, since the poten- 
tial decays slowly. Instead I choose to evolve the full 
domain of dependence of the initial data surface, hiding 
the effects of the boundary. 



IV. PARTICLE MOTION 

I use the same approach as described in paper II to 
evolve the particle's motion, i.e. I introduce the semi- 
latus rectum p, the eccentricity e and a fictitious angle 
X, not to be confused with the Faraday tensor component 
X defined in Eq. (|2.8p . such that 



r(r) = 



pM 



1 + ecosx{T) 
The evolution is then governed by 

dx (p — 2 — 2ecosx)(l -I- ecosx)^ 



(4.1) 



dt 



{Mp^) 



p — 6 — 2e cosx 
{p - 2 - 2e){p - 2 + 2e)' 



dip (p — 2 — 2e cosx)(l + e cosx)^ 
dt ~ p3/^My/{p - 2 - 2e)(p - 2 + 2e) 



(4.2) 



(4.3) 



I use the embedded Runge-Kutta-Fehlberg (4, 5) algo- 
rithm provided by the GNU Scientific Library routine 
gsl_odeiv_stepjrkf 45 and an adaptive step-size control 
to evolve the position of the particle forward in time. 



(5.1) 



ra=0 



where x = r* — Tq to the five points to the right of the 
particle's current position and extract ip and dr-tp as cq 
and Ci, respectively. In order to calculate ^^^^"'''"^ I fol- 
low [l3] and calculate '^^^^j^ ^^^'^ on the world line of the 
particle. Since this can be calculated using either the 
field values on the world line 



dHt,r*{t)) ^ 
dt 

ipit + h,r*{t + h)) - tp{t - h, r*{t - h)) 
2h 



(5.2) 



or as 



dHt,r*{t)) 
dt 



dip drQ 
dr*~dt 



where both and — tq/E are known, this allows 
me to find 



dip 
'dt 



d^{t,r*{t)) 
dt 



dip drl 
dr* dt 



(5.4) 



I repeat this procedure to the left of the particle. As a 
check for the extraction procedure, I compare the dif- 
ference between the right hand and left hand values 
[ip] — ?Aright — V'lcft with the analytically calculated jump 
conditions of appendix ID II Similarly I check whether 
the numerical solutions obtained for x and ^ directly are 
consistent with Eqs. p.l4p and ()2.13p . which give them 
in terms of derivatives of ip. 



VI. REGULARIZATION OF THE MODE SUM 



The regularization procedure operates on scalar spher- 
ical harmonic modes of the multipole coefficients F^l^^^^ 
of the Faraday tensor. As a first step I use the auxiliary 
fields Ip, X and ^ to reconstruct 



A 

dtv 



r,t 
I'm' 



(6.1a) 
(6.1b) 



and 



V. EXTRACTION OF FIELD DATA AT THE 
PARTICLE 

I use a straightforward one-sided extrapolation of field 
values to the right of the particle's position to extract 



a.t-^'™' = (6.1c) 

the combinations of the vector potential modes needed 
to obtain the even sector of a tensor spherical harmonic 
decomposition of the Faraday tensor. The auxiliary field 
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V and its derivatives provide the odd sector of the decom- 
position. 

Using the comples pseudo-Cartesian tetrad e^Q^, e"-|-) 
and e^gj introduced in paper I, I define tetrad compo- 
nents 

= Fa0e°'(t,)<^^(^) (6.2) 

of the Faraday tensor. 

I construct the spherical harmonic modes of 



using the couphng coefficients displayed in Eq. (jA4l 



F, 



£m,rct 



(.)(.) - E [qt)(.)(^'™'iM(<a"'-<;"' 

I'm' 



I'm' _ J^^' m' 



I calculate the multipole coefficients of ^(^^j^^) as 

Kit) = E<)(r;(^.'^o)F,™(|,v'o), (6.3) 

ra 

and regularize them as in Eq. (|1.5p . 



A 



B 



(6.4) 



(^-i)(^+|). 

I calculate the regularized self-force using F^^-^ — 
qF^^^^u^'^'^ . Finally I reconstruct the vector components 
of the self-force by from the tetrad components 



Ff = ^oF^,^, (6.5a) 
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FIG. 4: Convergence test of the numerical algorithm in the 
sourced case. I show differences between simulations using 
different step sizes of 16, 32 and 64 cells per M. Displayed 
are the rescaled differences 532-16 = £,{h — 1/32M) — ^(/i = 
1/16M) etc. of the field values at the position of the particle 
for a simulation with ^ = 6, m = 4 and p = 7, e = 0.3. I 
see that the convergence is approximately second-order. The 
curves are rescaled in such a way as to provide an estimate 
for the error of the highest resolution run compared to the 
real {h = 0) solution. 



A. Convergence tests 

Convergence tests are a straightforward way to test the 
implementation of a numerical scheme. I performed re- 
gression runs for my second-order convergent code using 
a non-zero charge q and an eccentric orbit. I extract the 
field at the position of the particle, and thus also test the 
implementation of the extraction algorithm described in 
section |Vl 1 choose the ^ = 6, m = 4 mode of the field 
generated by a particle on a mildly eccentric geodesic 
orbit with p = 7, e = 0.3. As shown in Fig. |4]the con- 
vergence is approximately of second order. In the region 
150 Af < t < 400 M the two curves lie on top of each 
other, as expected for a second-order convergent algo- 
rithm. In the region from 400 M to 450 M there is some 
difference between the two lines, caused by cell crossing 
effects similar to those discussed in paper II. 



VII. NUMERICAL TESTS 



In this section I present the tests I performed to vali- 
date my numerical evolution code. I performed the same 
set of tests as described in paper II. First, in order to 
check the second-order convergence rate of the code, I 
performed regression runs with increasing resolution. As 
a second test, I computed the regularized self- force for 
several different combinations of orbital elements p and e 
and checked that the multipole coefficients decay with £ 
as expected. This provided a very sensitive check on the 
overall implementation of the numerical scheme as well as 
the analytical calculations that lead to the regularization 
parameters. 



B. Discontinuity across the world line 

The singular source term on the right hand side of 
Eqs. (|2.15b|) - (|2.15f|) implies that the fields -0, x and £, 
are discontinuous across the world line. Since the jump 
conditions can be calculated analytically as done in ap- 
pendix ID 11 I can check whether the numerical results 
faithfully reproduce the expected behaviour. Using the 
methods described in section|V]l obtain one-sided extrap- 
olation for the field values and their spatial derivatives. 
For the highest resolution run used in the regression anal- 
ysis in section fVII Al I find that the numerical results for 
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FIG. 5: Multipole coefficients of the dimensionless Faraday 



■ImF, 



R 

(+)(-) 



for a particle on an eccentric 



tensor component 
orbit {p = 7.2, e = 0.5). The coefficients are extracted at 
t — 500 AI along the trajectory shown in Fig. [G] The plots 
show several stages of the regularization procedure, with a 
closer description of the curves to be found in the text. A 
uniform stepsize of h = 1/512 M was used. 



^ agree with the analytical calculation of the jump con- 
dition up to terms of the order of 10^^; two orders of 
magnitude smaller than the estimated numerical error of 
10^^. For dr*^ the situation is reversed, with the nu- 
merical error in the jump condition being about an order 
of magnitude larger than the numerical error in the field 
derivative itself. The accuracy of the numerical deriva- 
tives is therefore limited by the accuracy of the extraction 
scheme, resulting in about three significant figures for the 
set of parameters displayed in Fig. S) However the reg- 
ularization calculation is constructed in such a way that 
no derivatives of the fields need be obtained in order to 
calculate the self-force. I therefore feel that I can accept 
the reduced accuracy provided by the simple extraction 
scheme. 



C. High-^ behaviour of the multipole coefHcients 



Inspection of Eq. (|1.5p reveals that a plot of as 
a function of £ (for a fixed value of t) should display a 
linear growth in £ for large £. Removing the Aj^jjj^) term 
should produce a constant curve, removing the -B(;^)(iy) 
term (given that C(^)(y) = 0) should produce a curve 
that decays as and finally, removing the term 
should produce a curve that decays as £~'^. It is a pow- 
erful test of the overall implementation to check whether 
the numerical data behaves as expected. Fig. [5] plots the 
remainders as obtained from my numerical simulation, 
demonstrating the expected behaviour. It displays, on a 
logarithmic scale, the absolute value of ImF^^^^^^^^, the 
imaginary part of the F^.^^_.^ tetrad component of the 
Faraday tensor. The orbit is eccentric {p = 7.2, e = 0.5), 



and all components of the self-force require regulariza- 
tion. The first curve (in triangles) shows the unregular- 
ized multipole coefficients that increase linearly in £, as 
confirmed by fitting a straight line to the data. The sec- 
ond curve (in squares) shows partially regularized coeffi- 
cients, obtained after the removal of (^-|-l/2)^(^)(^); this 
clearly approaches a constant for large values of £. The 
curve made up of diamonds shows the behaviour after 
removal of B^^^-^^^^y, because C(^)(y) = 0, it decays as £~'^, 
a behaviour that is confirmed by a fit to the £>Z part of 
the curve. Finally, after removal of -D(^)(i,)/[(^— ^) (^+|)] 
the terms of the sum decrease in magnitude approxi- 
mately as £~^ when fitting to the data points £ > 1. 
This result depends slightly on the range of points used 
for the fit. I expect this to be due to the fact that I 
stop at £ — 15, which seems to be not large enough to 
show the asymptotic behaviour. Extending the range to 
very high values of £ proved to be very difficult, since the 
numerical code is only second order convergent, so that 
the numerical errors become dominant by the time the 
asymptotic behaviour begins to show. 

Each one of the last two curves would result in a con- 
verging sum, but the convergence is faster after subtract- 
ing the -D(^)(j/) terms. I thereby gain about one order of 
magnitude in the accuracy of the estimated sum. 

Figure[5]provides a sensitive test of the implementation 
of both the numerical and analytical parts of the calcu- 
lation. Small mistakes in either one will cause the differ- 
ence in Eq. (|1.5p to have a vastly different behaviour. 



D. Accuracy of the numerical method 

In this work I are less demanding with the numeri- 
cal accuracy then I were in paper II, where I describe a 
very high accuracy numerical code. Implementing suach 
a code is very tedious even for the scalar case, and much 
more so for the electromagnetic case treated here. There- 
fore I implement a simpler method that allows me to ac- 
cess the physics of the problem without being hindered 
by technical problems due to a complicated numerical 
method. 

An estimate for the truncation error arising from cut- 
ting short the summation in Eq. ()1.5|) at some -^max can 
be calculated by considering the behaviour of the remain- 
ing terms for large £. Detweiler et. al. [l^ showed that 
the remaining terms scale as £~'^ for large £. They find 
the functional form of the terms to be 



3/2 



{2£ - 3){2£ - 1){2£ + 3){2£ + 5) ' 



(7.1) 



where 7^3/2 ~ 36\/2. I fit a function of this form to the 
tail end of a plot of the multipole coefficients to find the 
coefficient E in Eq. (|7.ip . Extrapolating to £ —> 00 I find 
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error estimation mildly eccentric orbit 

relative truncation error in 2 x 10"'* 

relative discretization error in « 10"'' 

q ^ 



TABLE I: Estimated values for the various errors in the com- 
ponents of the self- force as described in the text. I show the 
truncation and discretization errors for the mildly eccentric 
orbit (p = 7.2, e = 0.5). The truncation error is calculated 
using a plot similar to the one shown in Fig. (5] The discretiza- 
tion error is estimated using a plot similar to that in Fig. H] 
for the ^ = 2, m = 2 mode. 

that the truncation error is 

oo 

e = [Eq. (dD] (7.2) 

{21 

max ~l~ 3)(2£iYLax H" l)(2^max l)(2-^inax 3) 

(7.3) 

where ^max is the value at vifhich I cut the summation 
short. 

A second source of error lies in the numerical calcula- 
tion of the retarded solution to the wave equation. This 
error depends on the step size h used to evolve the field 
forward in time. For a numerical scheme of a given con- 
vergence order, I can estimate this discretization error by 
extrapolating from simulations using different step sizes 
down to /i = 0. This is what was done in the graphs 
shown in Sec. IVII Al 

I display results for the mildly eccentric orbit shown in 
Fig. |6] with data extracted &t t — 500 M, that is at the 
instant shown in Fig. [S] At this moment, the multipole 
coefficients of Re(Fj^j) decay as expected, but e.g. the 

Im(i^j.^j) component decays faster with I for the range of 
modes < £ < 13 modes that were calculated. I choose 
an orbit of low eccentricity as high eccentricity causes 
the field values to be plagued by high frequency noise, as 
discussed in paper IL This makes it impossible to reliably 
estimate the discretization error for these orbits. 

Table H] lists typical values for the errors discussed 
above. 



VIII. SAMPLE RESULTS 

In this section I describe some results obtained from 
my numerical calculation. 

A. Mildly eccentric orbit 

I choose a particle on an eccentric orbit with p — 7.2, 
e = 0.5 which starts at r = pM/{l — e^), halfway be- 
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FIG. 6: Trajectory of a particle with p = 7.2, e — 0.5. The 
cross-hair indicates the point where the data for Fig. [S] was 
extracted. 
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FIG. 7: Regularized dimensionless self-force ^ij-Ft, -^-FV and 
on a particle on an eccentric orbit with p = 7.2, e = 0.5. 



tween periastron and apastron. The field is evolved for 
600 M with a uniform resolution of 512 grid points per 
M, both in the t and r* directions, for all values of i. 
Multipole coefficients for 1 < £ < 15 are calculated and 
used to reconstruct the regularized self-force Fa along 
the geodesic. Figure [7| shows the result of the calcula- 
tion. For the choice of parameters used to calculate the 
force shown in Fig. [71 the error bars corresponding to the 
truncation error Eq. (|7.2p (which are already much larger 
than than the discretization error) would be of the order 
of the line thickness and have not been drawn. 

Already for this small eccentricity, I see that the self- 
force is most important when the particle is closest to the 
black hole (ie. for 200 M <t<AOOM). The self-force 
acting on the particle is very small once the particle has 
moved away to r « 15 M. 
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FIG. 8: Trajectory of a particle on a zoom-whirl orbit with 
•p = 7.8001, e = 0.9. The cross-hairs indicate the positions 
where the data shown in Fig. 1101 and II II was extracted. 



B. Zoom- whirl orbit 



Particles on highly eccentric orbits are of most interest 
as sources of gravitational radiation. For nearly parabolic 
orbits with e < 1 and p 6 -f 2e, a particle revolves 
around the black hole a number of times, moving on 
a nearly circular trajectory close to the event horizon 
( "whirl phase" ) , before moving away from the black hole 
( "zoom phase" ) . During the whirl phase the particle is in 
the strong field region of the spacetime, emitting copious 
amounts of radiation. Figures |8] and [9] show the trajec- 
tory of a particle and the force on such an orbit with 
p = 7.8001, e = 0.9 calculated using a uniform step size 
oih = 1/256 throughout the range 1 < ^ < 15. Even 
more so than for the mildly eccentric orbit discussed in 
Sec. IVIII Al the self-force (and thus the amount of radi- 
ation produced) is much larger while the particle is close 
to the black hole than when it zooms out. The force 
graph is very similar to that obtained for the scalar self- 
force in paper II, however the overshooting behaviour at 
the onset and near the end of the whirl phase is not as 
pronounced. 

Since the rates of change in energy E and angular mo- 
mentum J of the trajectory are directly related to the 
self-force 



FIG. 9: Self-force acting on a particle. Shown is the dimen- 
sionless self- force ^Ft, and HFa, on a zoom- whirl 

orbit with p = 7.8001, e = 0.9. No error bars showing an 
estimate error are shown, as the errors shown are to small to 
show up on the graph. Notice that the self-force is essentially 
zero during the zoom phase 900 M <t < 1200 M and reaches 
a constant value very quickly after the particle enters into the 
whirl phase. 




E = -at, 



(8.1) 



FIG. 10: Multipole coefficients of — ReF(g) for a particle on 
a zoom-whirl orbit (p — 7.8001, e = 0.9). The coefficients 
are extracted at t = 525 A/ when the particle is deep within 
the whirl phase. Here r ~ and the behaviour of F^-^ ^ is 
very close to that for a circular orbit, requiring very little 
regularization. Red triangles are used for the unregularized 
multipole coefficients -F(o),£, squares, diamonds and disks are 
used for the partly regularized coefficients after the removal 
of the ^(0)) 5(0) S'lid D(o) terms respectively. 



it is easy to see that the self-force shown in Fig. [9] con- 
firms the expectation that the self-force decreases both 
the energy and angular momentum of the particle while 
radiation is emitted. 

In Fig. [TOl and Fig. [11] I show plots of F^q^ constructed 
from after the removal of the and 

terms. 



IX. EFFECTS OF THE CONSERVATIVE 
SELF-FORCE 

In this section only, I will use the subscript "0" to 
denote quantities evaluated on the unperturbed geodesic, 
and no subscript to denote quantities evaluated on the 
perturbed world-line. 
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FIG. 11: Multipole coefficients of — ReF^^ for a particle on 
a zoom-whirl orbit (p = 7.8001, e = 0.9). The coefficients are 
extracted at t = 1100 Af when the particle is far away from 
the black hole. As r is non-zero, all components of the self- 
force require regularization and I see that the dependence of 
the multipole coefficients on £ is as predicted by Eg. 11.51 After 
the removal of the regularization parameters Aj^jf^j, _B(^)(^), 
and -D(,i)(„) the remainder is proportional to £~'^ and £~'^ 
respectively. 

I follow the literature (see e.g. @) and define the dissi- 
pative part to be the half retarded minus half advanced 
force and the conservative part to be the half retarded 
plus half advanced force 



- { F' 



Fa) 



F 



adv\ 



(9.1) 
(9.2) 



The conservative force is the time reversal invariant part 
of the self-force. It does not affect the radiated energy or 
angular momentum fluxes E and J; it shifts the values 
of E and J away from their geodesic values, affecting the 
orbital motion and the phase of the emitted waves. 

To obtain expressions for E and J under the influence 
of the self- force, I employ the procedure described in [l6j . 
I begin by writing down the normalization condition for 
the four velocity 



1 



E^ 
f 



J2 



as well as the r-component of the geodesic equation 
F'' _ .. M .2 (r - 2M).P ME^ 



(r - 2M)r 



r4 (r-2A/)r' 
(9.4) 

where F"^ = qF^^u^ is the radial component of the self- 
force. Solving Eq. ([131) and (HH) I find 



, 9 (r - 2M)r F"" 
E^ =El-- ' 



3M m 



(9.5) 
(9.6) 



where 



. y (r-2M)rr (r - 2Mf 
Fl = r +— 771— + 71 TTTTZ' (9-7) 



V 



3M (r-3M)r' 
Mr2 



r - 3M 



3M 



(9.8) 



I stress that Eq and Jq are not the geodesic values for 
energy and angular momentum. They are of the correct 
form but are evaluated using the accelerated values for r, 
f and f (instead of the geodesic values rg, fg, etc.). 

For small perturbing force of order e I expand 
Eqs. (|9.5|) and (|9.6|) in terms of the perturbation strength 
and find 



E = Eo + AE Eq- e 



J = Jo + A J « Jo - e 



{r-2M)r F'' 
2{r - 3M)Eo~^ 

F-" 



(9.9) 



2(r - 3AJ) Jo m 



+ 0(£2), 



(9.10) 



where F^ is evaluated with the help of the unperturbed 
four velocity Uq = [Eo/ f,ro,0, Jq/tq]. The fractional 
changes AE/Eo and AJ/Jq are given by 

{r-2M)r F' 



AE/Eo 
AJ/Jo 



' 2{r - 3M)E^ m 

^4 pr 

' 2{r - 3M)J§~^ 



0{e' 



(9.11) 
(9.12) 



Once the perturbations in E and J are known, I cal- 
culate the change in the angular frequency 



'dt 



r - 2M J 
^ E' 



(9.13) 



For small perturbing force I expand in powers of the per- 
turbation strength 



ro - 2M Jo 
^ E'o 



l-e 



2{r- iM)Jl 
(r - 2M)r \ F ' 



2{r- 3M)E^ 



+ Ois'). (9.14) 



(9.3) The relative change Ai}/Qo is given by 



An/no 



{r-2M)r \F 



2{r - 3Af) J2 2{r - 3M)E^ J m 

+ 0(£2). (9.15) 



A. Circular orbits 

The effect of the conservative self-force is most clearly 
observed for circular orbits, where the unperturbed an- 
gular frequency fl as well as the shift due to the pertur- 
bation are constant in time. 
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FIG. 12: Fractional change Afi/fio induced by the presence 
of the conservative self-force. The effect of the self-force is to 
move the radius of the orbit outward, decreasing its angular 
frequency. 



For a particle in circular motion the self-force is con- 
stant in time and it turns out that the radial component 
is entirely conservative whereas the t and components 
are entirely dissipative. For circular orbits, the unper- 
turbed values of E and J are given by 

ro- 2M 

Eq = 



Jo 



^ro(ro - 3 A/) 



M 



ro -3Af' 

and substituting these into Eq. (I9.14p I find 



(ro - 3Af) [M 
2mAf 



(9.16) 
(9.17) 

(9.18) 



where the first term is just the angular frequency for an 
unperturbed geodesic at radius ro . The fractional change 
An/flo is then 



An 



(ro - 3Af )ro 
2mM 



(9.19) 



Similarly the fractional changes in E and J: AE/Eq and 
A J/ Jo are given by 



AE/Eo = -^Fr, 
2m 



AJ/Jo 



(ro - 2Ar)ro 



Fr 



(9.20) 
(9.21) 



2mM 

Figure [T^ shows the fractional change in J^q, E and J as 
a function of the orbit's radius ro. 



B. Eccentric orbits 

For eccentric orbits the self-force is no longer constant 
in time and I have to numerically calculate both the re- 
tarded and the advanced self-force in order to construct 
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FIG. 13: r component of the dimensionless self-force acting 
on a particle on a zoom-whirl orbit (p = 7.8001, e — 0.9) 
around a Schwarzschild black hole. Shown are the retarded 
(solid, red), advanced (dashed, green), conservative (dotted, 
blue) and dissipative (finely dotted, pink) force acting on the 
particle. 



the conservative self-force. I find the advanced force by 
running the simulation backwards in time. That is I start 
the evolution on the very last time slice and evolve back- 
wards in time until I reach the slice corresponding to 
t = 0. I reverse the boundary condition at the event 
horizon to be be outgoing radiation only (dt + dr* )V' = 
and adjust the outer boundary so as to simulate only the 
backwards domain of dependence of the initial slice. I do 
not change the trajectory of the particle. I do not change 
the regularization parameters, since they depend only on 
the local behaviour of the field and are insensitive to the 
boundary conditions far away. 



1. Conservative force on zoom-whirl orbits 

I calculate the conservative self-force on a zoom-whirl 
orbit with p = 7.8001, e = 0.9. Figs. [H and [H dis- 
play the breakdown of the self-force into retarded and ad- 
vanced, and conservative and dissipative parts for a par- 
ticle on a zoom- whirl orbit. In both plots the force is very 
weak when the particle is in the zoom phase t < 400 M 
or t > 800 M and nearly constant while the particle is in 
the whirl phase 400 M < t < 800 Af. Inspection of the 
behaviour of the r component reveals that it is almost ex- 
clusively conservative, with only a tiny dissipative effect 
when the particle enters or leaves the whirl phase. This 
result is consistent with the observation that the particle 
moves on a nearly circular trajectory while in the whirl 
phase, for which the radial component is precisely con- 
servative. Similarly the (f> component is almost entirely 
dissipative, with only a small conservative contribution 
when the particle enters or leaves the whirl phase, its 
maximum coinciding with that of f (not shown on the 
graph). 
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FIG. 14: if component of the dimensionless self-force acting 
on a particle on a zoom-whirl orbit (p = 7.8001, e = 0.9) 
around a Schwarzschild black hole. Shown are the retarded 
(solid, red), advanced (dashed, green), conservative (dotted, 
blue) and dissipative (finely dotted, pink) force acting on the 
particle. 



FIG. 16: r component of the dimensionless self-force acting 
on a particle on an orbit with p = 78, e — 0.9. Shown are 
the retarded and advanced forces as well as r. The vertical 
line at t ~ 2383 M marks the time of closest approach to the 
black hole. 
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FIG. 15: Relative change in Q, E, J for a particle on a zoom- 
whirl orbit due to the conservative electromagnetic self-force. 

I calculate the relative changes in J and under the 
influence of the self- force using Eqs. (ICT^l . (jO^ . ([051) . 
Fig. [15] displays the relative changes AE/Eq, AJ/Jq and 
Aft/riQ for a particle on a zoom whirl orbit p — 7.8001, 
e = 0.9. The change in E', J and il is strongest in the 
whirl phase when r 4.1M. It is consistent with the shift 
experienced by a particle on a circular orbit at 4.1M. 



C. Effects on the innermost stable orbit 

In the gravitational case, considerable work has been 
done to identify gauge invariant effects of the self- 
force [l3, [3 • The electromagnetic self- force is not sub- 
ject to the same ambiguity thus it can help shed light on 
the gravitational case as well by providing a clear distinc- 



tion between kinetic and dynamic effects. In this section 
I calculate the effect of the conservative self-foce on the 
location of the innermost stable circular orbit around a 
Schwarzschild black hole. Such a calculation was first 
performed for the scalar self- force by [l^ , where a highly 
accurate frequency domain numerical scheme was used. 
Recently 0, [3| have extended this calculation to gravity, 
using their time domain code to perform the intergration 
of the wave equation. Since the code presented in this 
paper is in the time domain as well, it is closest in spirit 

to \m]- 

X. RETARDATION OF THE SELF-FORCE 

For scalar perturbation in a weak gravitational field 
Poisson [23| showed that the self-force is delayed with re- 
spect to the particle motion by the light travel time from 
the particle to the central body and back to the particle 
again. In a spacetime where the central body is compact 
the treatment of [l^l is no longer directly applicable, but I 
still expect some retardation in the self-force when com- 
pared to the particle's motion. To study this effect, I 
calculate the self- force on an eccentric orbit with p — 78, 
e = 0.9; ten times larger than the zoom- whirl orbit dis- 
cussed earlier. The large orbit was chosen so as to be able 
to clearly see any possible retardation which might not 
be visible if the particle's orbit is deep within the strong 
field region close to the black hole. Figures [16] and [17] 
display plots of the r and (p components of the self-force 
acting on the particle close to periastron. Shown are the 
retarded and advanced forces as well as the particle's ra- 
dial velocity r. Without considering retardation I expect 
the self-force to be strongest when the particle is clos- 
est to the black hole, when r = 0, as evident in Fig. [T51 
Clearly for the r component displayed in Fig. [Tn]the re- 
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FIG. 17: (\> component of the dimensionless self-force acting 
on a particle on an orbit with p = 78, e = 0.9. Shown are 
the retarded and advanced forces as well as r. The vertical 
line at f ~ 2383 M marks the time of closest approach to the 
black hole. 
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FIG. 18: component of the retarded (solid red line) and 
(dashed green line) negative advanced self-forces acting on a 
particle with p = 78, e — 0.9. The forces have been shifted by 
At ~ 2[ro(t) — 3.0 M]. Also shown is the self-force calculated 
using the weak field expression Eg. 111. 1] (blue dotted line). 



tarded and advanced forces both peak at a time very 
close to the zero crossing of f, suggesting very little time 
delay in the r component of the self- force. In Fig. [TTlon 
the other hand the retarded and advanced <;i!)-component 
of the self-force peaks away from the time of closest ap- 
proach iniin- Inspection of the graph shows that the delay 
(advance) between the time of closest approach and the 
peak in the retarded (advanced) force is compatible with 
a delay of Ai„ii„ « 2(rmin - 3.0 M) « 74 Af. Using a de- 
lay of Ai « 2[Vo(i)-3.0M] and plotting F^^^t + At) and 
_^adv^^ — At) versus t both curves visually lie on top of 
each other and the maximum is located at imin as shown 
in Fig. [TH] below. This suggests that the self-force is in 
large parts due to radiation that travels into the strong 
field region close to the black hole and is scattered back 
to the particle. The time delay can then be loosely in- 
terpreted as the time it takes the signal to travel to the 
light ring around the black hole and back to the particle. 
This interpretation is loose for two reasons: First r* and 
not r is associated with the light travel time. Using r*, 
however, does not lead to a better overlap of the curves 
once a suitable constant offset chosen. Second, for the 
zoom- whirl orbit shown in Fig. 1^] the (shallow) maximum 
in the self-force is offset by only At sa 2[ro(i) - 1.0 M] 
which leads to a reasonable overlap of the two curves. 
Interestingly using r* instead of r yields a worse overlap. 
For very large orbits p = 780, e — 0.9 it is impossible to 
read off the small constant offset to the dominant 2ro(t) 
contribution. 



XI. WEAK FIELD LIMIT 

As a last application I use my code to compare the 
numerical self-force in the weak field region to the self- 



force calculated using the weak field expression 

fscli^Xc T'^ + ^rr- — — , g= J-r, (11.1) 

m r-^ 6 m dt 

of @, [llj . I calculate the self- force for a particle on an ec- 
centric orbits with e — 0.9 and p = 78 oi p — 780. Fig.[T51 
shows the retarded and (negative) advanced forces shifted 
by At « 2[ro(t) — 3.0 M] as well as the analytic force cal- 
culated using Eq. ()ll.ip . At this distance there are still 
some differences between the (shifted) retarded field and 
the weak field expression. One reason for this lies in the 
choice of a suitable r coordinate to correspond to the r 
coordinate in the weak field expression. In this work I 
use the areal Schwarzschild r, but the isotropic coordi- 
nate r = ^-^^ -\- even the tortoise r* could 
be used as well. Neither one yields a good agreement 
between the two curves. 

For p = 780 using a shift of At — 2ro (t) the agreement 
between numerical data and analytic expression is excel- 
lent as is evident in Fig. [1^1 At this distance r, f and r* 
are indistinguishable. 

XII. CONCLUSIONS 

I calculated the self-force acting a on an electromag- 
netic point charge in orbit around a Schwarzschild black 
hole. To do so I calculated the regularization parameters 
A,B, and D in sectionlBland implemented a second order 
accurate numerical scheme in section HIl 

I find the behaviour of the electromagnetic self-force 
to be similar but not identical to that of the scalar self- 
force. In both cases the self-force is strongest when the 
particle is closest to the black hole. Further, during the 
whirl phase of a zoom- whirl orbit with its nearly constant 
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Appendix A: Translation tables 



FIG. 19: (j> component of the retarded self-force acting on a 
particle on an orbit with p — 780, e — 0.9 close to periastron. 
Shown are the numerical (solid, red) and shifted analytical 
(dashed, green) forces. The agreement between numerical 
and analytical calculation is excellent, the discrepancy for t < 
7500 Ai is due to initial data contamination. 



radius, the self-force is very close to that of a particle in 
circular orbit at this radius. On the other hand, the 
overshooting effect upon entering the whirl phase which 
was observed in the scalar case is much weaker in the 
electromagnetic case. 

I calculated the effects of the conservative self-force on 
circular orbits, where it reduces the angular frequency 
and thus affects the phasing of the observed waves. I 
find this effect to be much stronger in the electromagnetic 
case than in the scalar case discussed in [l^ . In particu- 
lar during the nearly circular whirl phase of a zoom- whirl 

2 

orbit I find that decreases by ~ 0.06-;^. Due to the 

Mm 

2 

smallness of the ratio this change is tiny for one or- 
bit, however since it accumulates over the inspiral, its 
effect on the total phase shift during the full inspiral can 
be of order unity. This statement is not directly trans- 
ferable to the gravitational case since the radius ro of the 
orbit is not a gauge invariant quantity. Therefore I can- 
not distinguish between changes in due to effects of the 
self-force and due to gauge choices. To obtain a mean- 
ingful measure of the effect of the gravitational self-force 
I need to compare two gauge invariant quantities, e.g. 
and the gauge invariant ut of [22j . 

I investigated the retardation of the self-force with re- 
spect to the motion of the particle. I found that the 
retardation is very weak for the r component of the 
force and strong in the t and ip components, which are 
linked to radiated energy and angular momentum. In 
the later cases the retardation is compatible with a de- 
lay of At « 2(ro(i) — i?dciay)j where i?doiay is a constant 
depending on the particle's orbit. 



I require coupling coefficients to translate between the 
tensor harmonic modes of the Faraday tensor and the 
scalar harmonic modes of the tetrad components of the 
Faraday tensor. 

As a first step, I reconstruct the Faraday tensor modes 
from the numerical variables. For the even mode aux- 
iliary fields "0, X s-iid ? this reconstruction can be done 
algebraically while the odd sector requires a numerical 
differentiation of the numerical variable v. The recon- 
struction relations were already displayed in Eqs. (12. 9p - 
(|2.12p . which involves both the even and odd modes. 

In terms of the vector potential the Faraday ten- 
sor modes are reconstructed using the defining equation 
Eq. (12.11) . In this case, the reconstruction of the Faraday 
tensor reads 

Ftr = ^(4" - A^") (Ala) 
Ft A = J2 [("f" ~ ^*") ^a" + vT > (Alb) 

FrA = J2 [("f" ~ ^^") + > (Ale) 

l,rn 

Fe^ = Vim {X';} - X,^™). (Aid) 

Clearly both the expansion Eqs. (|2.9p - (I2.12p and the 
one in Eqs. (IAla[) - (IAld[) are of the same form and it is 
only necessary to obtain one set of translation coefficients 
to handle both the calculation using ip, x snd ^ in the 
main text and the one using the vector potential that will 
be presented in appendix [Cl 

The tetrad components are decomposed in 

terms of scalar spherical harmonics 

where each mode is given by 

^iJ^i^) - J n,K^)Yim dQ. (A3) 
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To obtain expressions for the coupling coefficients I sub- scalar case. 



stitute = ^a/3e"( )e'^. s into Eq. 



em 



\ a B 



+1? 



A': 



Ai 



(A4) 



which defines the coupling coefficients. It is often pos- 
sible to express these coupling coefficients in terms of 
linear combinations of the coupling coefficients derived 
in paper I for the scalar field. 

To simplify the notation of the coupling coefficients I 
use 



{l + ra){l + m+l) 
(2^+ l)(2£-h3) ' 

'(^ + m+ l)(£-m + 1) 
(2£+ l)(2£ + 3) ' 



(A5) 
(A6) 



as shorthands for recurring combinations of terms. With 
these the reusable scalar coupling coefficients are written 
as 

+ l'-"'+^^fh'l+l5ra'm-l (A7) 

r 

1^ + 2 



-7 



(A8) 



Ci^^){l'm'\lm) = -7' 



m'm+1 



y''^+^^-^5t,+i5^'m+i. (A9) 



E(^+){i'm'\£m) = — - + 1)(^ + H'^^'f^m'm-i, 

(AlO) 

i;(_)(^'m'|£m) = --v/(^ + m + 1)(^ - m)(5f.^5„,„+i. 

(All) 

In terms of the scalar coupling coefficients, the first 
coefficient for the expansion of i^(o)(-(-) is given by 



^ / y^'^'e^+^r^^dO (A12) 



1 

77 



and all other combinations of a, 5 and (/x) , {v) lead to van- 
ishing C^l^^^y Similarly for the remaining no n- vanishing 
coefficients for the -F'(o)(+) component 



D*o)(^)(fm'|M = C(+)(/m'|M/V7, (A13) 
i?[o)(+)(^'m'|£m) = E^+^{£'m'\£m)/^. (A14) 

The coupling coefficients for contain both 

even and odd modes. The first non-vanishing one is 
DJ^_^-^^_-^{£'m'\£m), which is given by 

m'Km) = v/77''"™+'C^(-)(^'™'K + 1, m - 1) 

- Vll'-''"'C(^^)i£'m'\£~ l,m- 1) 

+ /77''"+'q+) ii'm'\£ + 1, m + 1) 

_y7^^-i,-™C(+)(fm'K-l,m + l), 

(A15) 

while the coefficients coupling to odd modes are 

El+)^_){£'m'\£m) = ^/fY''"'+^E^_^{£'m'\£+l,m-l) 
-/77'^'''"£^(-)(^'™'K-l,m-l) 

-^/fl'-'^-"'E^+){£'m'\£~l,m + l), 

(A16) 

i;(+)(_)(/mVm) = -2i/r\£ + l)i£ + 2)e'"'5t,t+i5m'm 



2i/r'{£-l)£e''-'^"'5t't-i5. 



(A17) 



Appendix B: Regularization parameters 



Similarly it proves useful to define lower order coupling 
coefficients for the odd sector, which is absent in the 



I mimic the treatment in paper I and start from the co- 
variant expression for the singular vector potential tensor 
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in Eq. (464) of [23^ 

+ TrUaP'-ISu'^ + 7rUal3'--i'u'^ U' V/3U+ UaP"-l3u'^ 

ZT ZT -^^adv 

1 // 1 

ypv--q / V/3V^a^(a;,z(T))u^(T)dT. 

(Bl) 

I have introduced a large number of symbols, x is the 
point where the field is evaluated, x' and x" are the re- 
tarded and advanced points of x on the world line z{t). 
They are connected to x with unique future-directed and 
past-directed null geodesies, respectively. u(x) and v{x) 
are the retarded and advanced time functions such that 
x' — z{t = u), x" = z[t — v). u°' and are the 
four velocity at x' and x" respectively. Further I define 
Synge's world function <t{x, x) which is numerically equal 
to half the squared geodesic distance between two points 
X and X. Using its gradient = "S/ ci'^{x,x), I define 
r = u'^ <7a'{x,x') and radv — — aa"{x,x"), the afhne 
parameter distances of x away from the world line along 
its future/past light cone. The potentials U and V ap- 
pearing in Eq. (jBip are the direct and tail parts of the 
retarded Green function Gi^p{x,x) associated with the 
wave operator. 

Prom the definition of r, radv, u, and v it follows that 
(see Section 3.3.3 of [13) 

V„u = -(T„(x,x')/r, (B2) 
VaV = aa{x,x")/r^dv, (B3) 

Va?' = —(Ja'P'U" U*^ VttU + (Ta'aU"' , (B4) 

V^radv = -a'Q"/3"it" u'^ VqM - o-q"qu" , (B5) 

which are valid for geodesic motion. 

The potentials UaiS', Uap" are determined by Eq. (322) 

of m 

C//' =/>i/2(a;,a:'), (B6) 
C//"=/>i/2(x,x"), (B7) 

where g'^^ is the parallel propagator from x'^ to x^ and 
A = det(— ^crf^,) is the van Vleck determinant. Of 
the potentials Vap' and Vap" appearing in Eq. (|B1[) I 
only need to know the scaling behaviour following from 
Eq. (320) of dl: 

14/3' = O(e'), (B8) 
V^p„ = 0{e^), (B9) 
V^y^^ = 0{e). (BlO) 

These expressions are valid in vacuum spacetimes where 
the Ricci tensor vanishes. 



Again mirroring the calculation in paper I, I introduce 
the arbitrary point x = z(f ) on the world line and expand 
the quantities in Eq. ()B1[) in terms of a Taylor expansion 
around x. I introduce the convenient quantities 



r = aa{x, x)u°' , (Bll) 

s2 = (g"^ -t-u"u^)cra(a;,x)cr^(a;,5), (B12) 

together with the time differences 

A+=w-f, A_ = u-f (B13) 



from the advanced (retarded) point to the reference point 

X. 

I also use the expansion of the derivatives of the par- 
allel propagator around the point x 

(B14) 

-fO(e^), (B15) 

as well as an expansion for the second derivative of 
Synge's world function 

+ ^Ra^Bh^^^^'^' + 0(£4), (B16) 

and the van Vleck determinant 

Ai/2 ^ l + 0(e4), (B17) 

which I calculate using the methods described in 
Sec. (2.4.2) of 

I make use of the fact that the bi-tensors 

C/a(r) = [/,X- (B18) 
Ua.p{T) = Ua,^,pu>',&nd (B19) 

Uc.{t) = C/a^;,u'^u'' (B20) 

appearing in Eq. (IB1[) do not bear a free index on the 
world line, making them scalars on the world line. With 
w being either m or w and A = w — f = A^ I expand 
these as 

Uo.{w) = + f7„A + ij7,A2 + ^{/^=^)a3 + 0{e^), 

(B21) 

Uo.0{w) = Uo.p + f/a,3A + ^Uc^pA^ + 0{e^), and (B22) 

Ua{w) - f7„ + UaA + ic/^^^A^ + 0(e3), (b23) 

where it is understood that the coefficient functions are 
evaluated at t = f. 
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Repeatedly taking derivatives of Eq. (jB7l) and con- 
tracting with u'^ I find for the first set of coefficients 



(B24) 



I, +0(£^), (B25) 



Uc.^ p\Rauua\u+0{e^), (B26) 
[/^3) =0 + O(£), (B27) 
where I have introduced the notation Raima = 



Rap^gu^u'a^ and R, 



aUU(T\(T 



R^^^s.^.u^u^a'a^-; I will 
use this notation and its natural extension to higher 
derivatives and different combinations of u'' and a'^ fre- 
quently below. 



Similarly I find for the second set 

3 



UaP = g^-y^ ( ^Rs.uPa - ^R^uf,a\a ] + 0{e'), (B28) 



UaP — 9"a9^ a ( T^Raufiu + ^RauPalu ~ nRauBu 



1 



1 



^ V2 



6 



U.l3u\(T 



(B29) 



Uap = ^9\9%Rs.uPu\u + 0{e). (B30) 

Note that the third set does not involve new coefficients, 
but only those already calculated for Ua- 

Finally I copy expressions for A±, r, radv, u, v and 
their gradients from paper I 



A± = (f ± S) T Ruaua T [(f ± S)i?„^„^|„ - Ruaua\a\ + O(e^), 

-9 9 



r = s 



^uau(j\u 



f adv = S - 



2 - 

-^i?n^ua-^^ [[r- s){f + 2s)R 

u 

'uiTua 24s~ '^'^^^ ~ '^s)RuaucT\u ~ ~ s)-RMcruCT|CT] + 0(£^), 



J.2 _ g2 



6s 

CTa + (f - s)ua 



i i 

— (^r SjltatTua 77 

6 3 

^)-^Qiicrii|a" 



r — s 



((r - s)(r -I- 2s)i?„^„^|„ - (r -|- s)i? ) CTS 



24s2 



"1 l._ .2 r^-s^ (f + s)2(f-2s) 



(B31) 
(B32) 
(B33) 



(B34) 



(B35) 



19 



V^r = g"^ I [o-Q + rua] 

s ' 



1_ 2n f^+S 

6 3 6s"^ 



fo2 



1 



— (r - s)2(r + 2s)Rs,uau\u + {{r - s){r^ + fs + 4:S^)Ruaua\u - {r^ + s^)Ruaua\a) 0-a 



12 

f — s 



+ (('^ - s)ir^ + 2rs + 3s^)Ruaua\u - rir + s)i?„^„ak) 



0(6^) , 



Var-adv = g"a \ i'^a + ''"a] + 



1_ 1 2 2^ f'^+s'^ 

H 

^6 3 bs'^ 

-—fR- I +iff2-s2)i?- I +—(f^-s^)R- I 



1 

12' 
f + s 



1 



^To + ^ '2.s)Rauau\u + (l'' + s){r^ ~rs + ^S^)Rucrua\u - [r^ + S^)Ruaua\a) 



■U(7ua'\u 



{t s')Ruau(T\(T^ 

I 



(B36) 



(B37) 



After substituting Eqs. (jB2p - (jB37P into Eq. (|Bip (all pression for the covariant expansion of 
of them) and sorting out the orders I find the final ex- 



-9 2 



T / V . ^ 



3f2 - s2 



3s2 
f 



r{r'^ — s^) 
2^5 



r 

RuauijUaUp + 'l^R-aufiu 



j,2 _ g2 



^2 ^2 \ ^2 ^2 



1 

'3s 



_ p 

2 au/3(T|(T 



6s 



fo3 



r(r2 - 3s2) 



r 



6s 



-R. 



aubcrlu 



^^3 ^aiibMlii 



j Rauua\u ~l" 12s'^ ^OiRfiuua\u -^Q;u/3u|dr 

3^2 — f{f^ — s^) 



24s5 

7~(f2 — S^) 



g^5 -^■^uaua\u j 13 
(f2 - s2)2 



8s5 



-0{e^), (B38) 



r 



where terms in square brackets are of the same power in 
e. 

I copy the results for the coordinate expansion of 
aa{x,x) and g°'p{x,x) from Eqs. (3.16) - (3.19) and 
Eqs. (3.30) - (3.33) of paper I. I use 



-aaix, x) = gapW^ + Aap-fW^w'^ + AalS^SW^w'^l 



Aa — 1 /pa I T^a ^ 



(B39) 
(B40) 

(B41) 



24 (^"''^^ 



I pa pP 



lOpa p/i I pct pA" pi' 



as well as 



g\{x,x)^6^^^ + B\y +B^^^ 



Pi 



W' W 



a/37 2 



2 V 



T^M pi' 

7 ' 01/ a7 



(B42) 



(B43) 
(B44) 

(B45) 
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Pl^ ab,v ~ ^ fiv^ qAJ- 75 

+r'^,^r^„^r\, + 2r''^,r%^r\,) , (B46) 

where = — a;" is the coordmate distance between 
X and X. Together with Eq. (jB38P these equations form 
an expansion in of the singular part of the gradient of 

I 



the vector potential around a point x near the world line 
of the particle. I finally calculate the tetrad components 
of the singular Faraday tensor as 

From this point on I proceed exactly as described in Sec- 
tion V of paper I using Maple and GRTensorII to per- 
form the calculations. I find, after an extremely tedious 
calculation, 



^(o)(+) = sign(A) 



(o)(+) 



Trfaro 



D 



(o)(+) 



iEr^i-Urp^ + J* + r^)fg (-rof + 2ro + 7r^f - Ur^)Er\ 



+ 1 



{SMrofJ^ - SrgMJ" + SSJVjjf - 2rg - WMr^J^ + SJV^f + 5U^r^f + 20r^f + 5r§f - 6r^)E 



iE4{lJ^-rl)rl [2-f)TlEfl 
(AMrgfJ^ + 20.Prf,f + 8Mrp^ + 14r^Mf - 2rg + 12Afrg + 4f)E 

i{2MJ^ - Mlrp^ - 2 jVgf - 20Mrp^ - 2r5f + 2r5 - lMIrl)Era 



4rj]a5fj7r 



^(+)(-) = sign(A) 



2iEJ 



B 



(+)(-) 



-2i< 



jrl ~ J^)fl ^ -J^ro f + 2ro,P + 2rl - 2rlp 



r'or'o 2(1 - f) 



^(+)(-) =-2«< 



43 jV;5fA/ - 7jV;^f - 27J^Mrl^ - llAfrgf - r-jf + r|j - 2rgA/)ro/^f^ / (4r5a-Vf7rj 
(sA/f - 8 J^Af + 30Af rgf - 2 jV^ + 10 jVj^fAf - 24 jV^^Af + 3 jV^f 

28 J^Az/r^ + jVJf - 28 J^Af r|;f - 20A/r^f - 12r^Af) j (srla^Jir 
8f7rJa7 

4Af f J6 - 16 J^A/rg + Arl^MJ^ ~ 18 J^r^JfAf - 28 J^rfjA// - J^fr^ - 12r§A4" - 20r3fAf 
(2MfJ^ ~ ^rl]M3^ + J'^rl - 2 jV^Af - 5 J^fr^ - MJ^r^fM - 2r^A/ -t- 



llr 



;fA.f-fr^)foy(4ro^/2aVf7r 



(B48) 
(B49) 



(B50) 
(B51) 
(B52) 



(B53) 
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where f=.y^^, a' = rl + j\ 

Here, the rescaled elhptic integrals 8 and JC are defined 
by 

f = - (1 - fcsin^ ^fl^ dV' = F f -; 1; fcj , 

(B54) 

and 

/C=-y (l-A;sin2^)-i/2d^ = Ff-,-;l;fcj , 

(B55) 

in which i^(a, &; c; a;) are the hypergeometric functions 
and fc EE J7(r2 + J2). 

Appendix C: Vector potential calculation 

In this section I describe a variant of the numerical 
calculation discussed in the main part of the paper that 
uses the vector potential instead of the Faraday tensor. 
To this end I decompose the vector potential and the 
sources in terms of vectorial spherical harmonics 

A,(i,r*,0,0) = iAf^(i,r*)y,„(0,0), (Cla) 

ja(t, r\Q, 4>) = fj^it, r*)Yem{e, 0) for a^t,r*, 

(Clb) 



I 

AA{t,r* ,0,4,) ^ ve,n{t,r*)Z'r{0,c^) 

+ ve^{t,r*)X'r{e,4,), (Clc) 

+ (t,r*)Xl"(0,0) for A = 0,0, 

(Cld) 



and substitute this into the Maxwell equations for the 
vector potential in the Lorenz gauge g^^Aa-p = 0: 



g^'' A^,^, ^ R\A^ = -A7:j^, (C2) 



where Rap is the spacetime's Ricci tensor, which vanishes 
in Schwarzschild spacetime. Substituting Eq. (IC1|) into 
Eq. (|C2|) I arrive at two decoupled sets of equations for 
the even (A^™, vem) and odd {vim) modes 



2M fdAil" 

~ \ ^ 

dr* 



dr 



*2 



dr* r 
2M /dAf"^ 



dr* 

-(v. 



Vvir, 



2^1 AiT 



fVvim = 



-Anrff^T, 



2^Ai^ = -Anfj, 



even 
'm ' 



dt^ 



(C3) 

(C4) 

(C5) 
(C6) 



where V and j^™ is defined as in Eqs. (|2.15bp and (|2.4p 
in the main text. 



1. Numerical method 

I discretize the set of reduced equations Eqs. (jC3| - 
()C6P using Lousto's method as described in section ITTl of 
the main text. Since the source terms on the right hand 
side are less singular for the vector potential than they 
are for the Faraday tensor, I do not have to distinguish 
between sourced and vacuum cells in the integral over 
the potential terms. 



Terms containing first derivatives ^ , , where now 
and in the remainder of the appendix ij) stands for any of 
were not treated in [lo|, but, for 



Aim 



A. 



or V , 



generic vacuum cells, can be handled in a straightforward 
manner 



dudvV{r)^ = 2h{^j3 - ij2)Vo + 0{h''), (C7) 
cell c/r 

dudvV{r)p^ = 2h{^i-^ji)Vo + 0{h^). (C8) 

cell C'''* 



This fails for cells traversed by the particle, since the field 
is only continuous across the world line but not differen- 
tiable. For these cells I take recourse to Lousto's original 
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algorithm, which has to deal with a similar issue, and use 



4. Initial values and boundary conditions 



dudvV{r)^ ^VoJ2AdtiP, + Oih''), (C9) 

i 

dudvV{r)^ = Vo + (CIO) 



where Ai,. . . ,^4 are the subareas indicated in Fig. [3] and 
dt4'i, ■ ■ ■ , dtip4, dr'ipii ■ • ■ : dr*ip4: are zeroth order accu- 
rate approximations to the derivatives in the subareas. I 
calculate these using grid points outside of the cell on the 
same side of the world line as the corresponding subarea, 
e.g. 

d..^, = 'f'^''^*-'^-J^'''*~''Koih). (Cll) 



2. Gauge condition 



I handle the problem of initial data and boundary con- 
ditions the same way as in the main text, that is I ar- 
bitrarily choose the fields to vanish on the characteristic 
slices u = uo and v = vq 



Aa{u = Uo) = Aa{v = Vq) = 0, 



(CM) 



thereby adding a certain amount of spurious waves to 
the solution which show up as an initial burst. Gauge 
violations in this initial data are damped out along with 
those arising during the evolution. 

I implement ingoing wave boundary conditions near 
the event horizon and choose a numerical domain that 
covers the full domain of dependence of the initial data 
near the outer boundary. 



In contrast to the scalar field, the electromagnetic vec- 
tor potential has to satisfy a gauge condition 



g"^A^,f3 = 0. 



(C12) 



Analytically the gauge condition is preserved by the evo- 
lution equations, so that it is sufficient to impose it on 
the initial data. Numerically, however, small violations of 
the gauge condition due to the numerical approximation 
can be amplified exponentially and come to dominate 
the numerical data. To handle this situation I introduce 
a gauge damping scheme as described in [23, HI] . That 
is I add a term of the form 



5. Extraction of the field data at the particle 

In order to extract the value of the fields and their 
first derivatives at the position of the particle, I use a 
variant of the extraction scheme described in paper II. I 
introduce a piecewise polynomial 



p{x) 



Co -I- Cix + if a: < 



Co + CiX + 



if a; > 



(C15) 



1 



dAt 



1 



dAr 



-2M dt 



1 , 



r-2M dr 



(C13) 



to the t components of the evolution equations Eqs. (jC2l) . 
which dampens out violations of the gauge condition. 
This choice proved to be numerically stable for the radia- 
tive {£ > 0) modes but unstable for the monopole {£ = 0) 
mode. 



in a; = r* — Tq on the current slice. Its coefficients to the 
left and right of the world line are linked by jump condi- 
tions c„ = c'ji + [d",ip] listed in Appendix lD 21 Fitting this 
polynomial to the three grid points closest to the parti- 



cle, I extract approximations for ^p{to,rQ) 



and £%ri2 

or 

which are just the coefficients cq, ci respectively. Once 
I have obtained these, I proceed as in section |V] of the 
main part of the paper following [T3| to obtain values for 

dt 



3. Monopole mode 

The monopole moment of an electromagnetic field is 
non-radiative. This makes its behaviour sufficiently dif- 
ferent from that of the radiative {£ > 0) modes that the 
approach outlined earlier fails for £ — 0. In this case 
Eq. (|C2p reduces to a set of coupled equations for 
only. Rather than solving the system of equations di- 
rectly for and A^ I use the analytical result for 
the Ftr component of the Faraday tensor derived in sec- 
tion III CI in the main part of the paper. This jjroves to 
be sufficient to reconstruct the combination A^.^ — A^ \ 
appearing in Eq. (jA4p . 



6. Results 

Using the vector potential code described above I can 
reproduce the results obtained from the Faraday tensor 
method discussed in the main paper. The differences are 
small, typically of the order of 10~^% of the field values 
as shown in Fig. [20l I expect the Faraday tensor code 
to yield more accurate results since the costly numerical 
differentiation that is necessary in the vector potential 
calculation is absent. Nevertheless I can reproduce e.g. 
the correct decay behaviour of the multipole coefficients 
for a zoom- whirl orbit as shown in Fig. 1211 
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FIG. 20: Differences between calculated using the vector 
potential and calculated using the Faraday tensor method for 
£ = 2, m = 2 mode of field for the zoom- whirl orbit shown 
in Fig. [8l Displayed are the difference and the actual field. 
The stepsizes were h = 1.0416 x 10"^ M and h = 1/512 M 
for the vector potential calculation and the Faraday tensor 
calculation respectively. 
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FIG. 21: Multipole coefficients of — ReF^Qj for a particle on 
a zoom-whirl orbit (p = 7.8001, e — 0.9), calculated using a 
stepsize of /i = 0.125M for the £ — 1 modes and increasing 
the resolution linearly with £ for £ > 1. The coefficients are 
extracted a,t t = 1100 M when the particle is deep within 
the zoom phase. Red triangles are used for the unregularized 
multipole coefficients f"(o),«) squares, diamonds and disks are 
used for the partly regularized coefficients after the removal 
of the ^(0), S(o) and _D(o) terms respectively. 



Appendix D: Jump conditions 



1. Faraday tensor calculation 



Since the source term in Eqs. (j2.15bp - (j2.15f|) contains 
a term proportional to S'{r* — Tq), the field is discontin- 
uous across the world line of the particle. I use 



6^0+ 

-d^d^,^P{to,r*~e)] (Dl) 

to denote the jump in d^d^ip across the world line. I 
only calculate jump conditions in the r* direction up to 
which I find by substituting the ansatz 



(D2) 



into Eqs. (|2.15bp - (j2.15f|) and its t and r* derivatives. 
Demanding in each step that the singularity structure 
on the left hand side matches that of the sources (and 
their derivatives) on the right hand side yields the jump 
conditions 



and 



/oP.r5)2-l]' 



G,i 



(D3) 



/o[(5,r*)2-l]3 

where ip stands for either one of V', X, or ^. 

2. Vector potential calculation 



(D4) 



Since the source term in Eq. (|C2[) is singular, the field 
is only continuous across the world line of the particle, 
but not smooth. I use 

[^r^^V-] = hm [ara-V(io,rS +£) 

e-^0+ 

-d:d^Mto,r*o-e)] (D5) 

to denote the jump in dj^d^ip across the world line. For 
my purposes I only need the jump conditions in the r* 
direction up to [9^* ?/'] i which I find by substituting the 
ansatz 



(D6) 
(D7) 
(D8) 



into Eqs. (jC3l) - (jC6[) and its t and r* derivatives. De- 
manding in each step that the singularity structure on 
the left hand side matches that of the sources (and their 



v''"'{t,r*) ^ v''^{t,r*)0{r* - r*) 

-)-<(t,r*)0(r*-r*), 
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derivatives) on the right hand side yields the jump con- 
ditions 



[dr. A': 



w 



E^ 



• 2 *^cvcn/odd; 



(D9) 
(DIO) 

(Dll) 



2ME'^ 



2(R2 _ ^2\2 



rl{E 



(3rl + E^)EW, . 



2ME^r 



rliE^-PoY 



{E^ 
Sb — fo 



r'5'a 



2E^ro 
[E^-flY' 

(D12) 



[5; 



= -/o 



{Zrl + E^)E^r^ 
2E^ro 



Sr 



/( 



So: 



n/odd 



i/odd? 



where a^b G {t. r*}, a ^ w G {i;, v}. 



(D13) 
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